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ABSTRACT 


A prototype All-Reflection Michelson Interferometer (AMI) was tested using two 
different diffraction gratings and a variety of light sources. The original design used a focusing 
lens, a pinhole aperture and an off-axis parabolic mirror for collimation, a plane diffraction 
grating and two plane mirrors to divide and recombine the light and a CCD camera to record the 
interference pattern. Interference patterns were recorded and analyzed in the Interactive Data 
Language (IDL) using Fourier transform techniques. The design was then modified to 
approximate an actual instrument for remote observations of the atomic oxygen triplet emission. 
The results of numerous experiments using both the original and modified instruments verified 
the feasibility of using the AMI as a compact, lightweight, high resolution instrument for use on 
sounding rocket or satellite platforms. 
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I. INTRODUCTION 


Knowledge of the chemistry and composition of the earth's upper atmosphere is of vital 
importance to many academic and government activities. Specifically, an understanding of the 
electron density profile of the ionosphere has direct application to military over-the-horizon 
radar systems and high frequency radio communications. Understanding the thermosphere is also 
important to the scientific community in areas ranging from accurate prediction of weather 
patterns to the study of terrestrial planets and comets. 

Atomic oxygen is a major constituent of the upper atmosphere of the earth and, because 
it is highly reactive, plays a central role in the photochemistry of the ionosphere. It is formed 
naturally in the thermosphere by the photo-dissociation and electron-impact dissociation of 
molecular oxygen. Employing a variety of methods, the density of atomic oxygen in the 
thermosphere has been measured on numerous occasions, but the results differ widely and the 
absolute concentration of atomic oxygen has proven notoriously difficult to establish. 

A relatively new technique for measuring atomic oxygen in the thermosphere under 
research at the Naval Postgraduate School (NPS) involves using the solar occultation of the OI 
1304-A emission to infer the terrestrial atomic oxygen density. This technique is unique because 
it is insensitive to many of the sources of uncertainty that have affected previous measurements. 
The OI 1304-A triplet is a prominent emission in both the solar spectrum and terrestrial airglow. 
Any observation of the solar OI 1304-A line profile from within the atmosphere will be absorbed 
to some degree by the intervening column of oxygen and the extent of this attenuation can be 
used to determine the column density of the oxygen. 



For space flight observations, high resolution measurements are made almost exclusively 
with Michelson or Fabry-Perot interferometers rather than grating spectrometers. Emissions at 
wavelengths of 1304-A fall into the far ultraviolet (FUV) portion of the spectrum. The number of 
window materials in the FUV is very limited and, as a result, Fabry-Perot and standard 
Michelson interferometers are extremely difficult to build. Therefore, instruments utilizing all¬ 
reflection interferometers are needed for this portion of the spectrum. 

The NPS Physics Department has been actively involved in ionospheric research for 
several years. They have conducted three successful sounding rocket experiments using an 
instrument called the Middle Ultraviolet Spectrograph (MUSTANG). MUSTANG observes 
emissions in the range from 1800-A to 3400-A with a resolution of approximately 10-A. The 
successor to MUSTANG is the higher resolution Ionospheric Spectroscopy And Atmospheric 
Composition (ISAAC) instrument, scheduled for launch in 1995 aboard a low earth orbit 
satellite. ISAAC is expected to have a resolution of 2.5-A. 

During MUSTANG development it became evident that an instrument providing even 
higher wavelength resolution, on the order of a few mA, was necessary to permit the 
investigation of radiation transport through the examination of optically thick lines. Nichols 
(1990) designed and built a prototype interferometer using all-reflective optics and a spherical 
grating. The design was validated using a green He-Ne laser in the laboratory (Cleary et al., 
1992). Wallace (1992) found that a similar instrument based on a plane diffraction grating is 
more suitable because it produces linear interference fringes which are more amenable to 
analysis using Fourier transform techniques. He designed and built the All-Reflection Michelson 
Interferometer (AMI). Validation of the design using the DART ray tracing program was 



performed by Carlson (1993) and developmental testing was begun by Risley (1993). This thesis 
is a continuation of that work. 

A. THESIS OBJECTIVES 

The main objective of this thesis is to continue developmental testing of the AMI using 
visible light sources in the 5000 to 6000-A range and incoherent light sources with wavelengths 
in the 2000 to 3000-A range. The ultimate goal is to demonstrate the feasibility of the AMI for 
analyzing the solar occultation of the OI 1304-A triplet emission to determine the density of 
atomic oxygen in the thermosphere so that development of a practical instrument for space flight 
measurements can begin. 

B. THESIS OUTLINE 

This thesis is divided into five sections. The first section is the introduction of the thesis. 
Section II provides background information on the composition of the atmosphere and the 
dynamic processes that produce the emission profiles of interest. The theory of Section III 
includes a brief review of interference and diffraction, an introduction to the AMI and solar 
occultation. Section IV describes the experiments conducted and data collected from the various 
instruments, and Section V presents the results of the research and provides recommendations 
for further study. Appendices are also included which contain IDL procedures used to store and 
to analyze the data and EDC files used to generate some of the figures. 



H. BACKGROUND 


This section provides general information relating to the physical processes and 
dynamics which control the composition of the earth’s atmosphere. It also provides the 
motivation which led to the development of a high resolution FUV interferometer. 

A. ATMOSPHERIC COMPOSITION OF THE EARTH 
1. Temperature and Composition Profiles 

The earth’s atmosphere can be described as a series of stacked horizontal layers 


based on kinetic temperature or chemical composition. The temperature and composition 
distributions of the atmosphere are illustrated in Figure 1. 





Temperature Profile of the Atmosphere 





The atmosphere is divided into four different regions or “spheres” based 
primarily on the variation of temperature with altitude above the earth’s surface. The different 
regions in order of ascending altitude are: the troposphere, the stratosphere, the mesosphere and 
the thermosphere. Each of these spheres is separated by a transition region or “pause” where the 
derivative of temperature with respect to altitude goes to zero. The three transition regions 
separating the four characteristic layers are the tropopause, the stratopause and the mesopause. 

The troposphere (0-10 km) is the layer of the atmosphere where clouds and 
weather originate. In this region, the dominant physical process is Eddy mixing and the surface 
of the earth is the primary source of warming. As a result, the temperature decreases at a fairly 
uniform rate (commonly referred to as the lapse rate) of -6.5 K/km and reaches a minimum of 
approximately 210 K at the tropopause (10 km). 

The next region of the atmosphere is the stratosphere (10-50 km), where the 
temperature increases with increasing altitude due to the absorption of solar ultraviolet radiation. 
This effect is especially true in the ozone layer (30 km) where a concentration of a few parts per 
million of 0 3 leads to very pronounced absorption of wavelengths in the 2000-3000-A range and 
subsequent heating of the atmospheric gases. The positive temperature gradient observed in the 
stratosphere terminates at the stratopause (50 km). 

The third layer is called the mesosphere (50-80 km) which is characterized by a 
negative temperature gradient. The temperature in this region is determined mostly by 
convection, the warmer stratosphere now acting in much the same way as the planetary surface 
acts with respect to the troposphere. The cooling trend observed in the mesosphere ceases at the 
mesopause (80 km), the coolest region of the atmosphere, where the temperature is 
approximately 180 K. 



The Final region of the atmosphere is the thermosphere (above 80 km), which 
displays a positive temperature gradient due primarily to the absorption of electromagnetic 
radiation by molecular oxygen in the range of X< 2000-A. Since the heating now occurs from 
above, there is no convection, and the temperature increases until it reaches a maximum of 
approximately 1600 K at altitudes of 300-400 km. Above 400 km the atmosphere is essentially 
isothermal. As Figure 1 shows, the high temperatures found in the upper reaches of the 
atmosphere display considerable variation (900-1800 K) depending on the amount of solar 
activity and time-of-day. These variations are a result of the low heat capacity of the upper 
atmosphere. Fluctuations in the amount of UV radiation are rapidly reflected in temperature 
changes. In the thermosphere, interaction between solar radiation and atmospheric gases form 
the region of ionized particles known as the ionosphere. This region is of primary interest for 
applications discussed in this thesis. 

If the atmosphere is divided according to composition there are three 
characteristic spheres: the homosphere, the heterosphere and the exosphere. In the homosphere 
(0-100 km), the composition of the atmosphere is essentially uniform in the sense that the mean 
molecular weight remains constant. The primary transport mechanism in this region is eddy 
diffusion, a turbulent mixing which keeps the relative concentrations of component gases 
effectively constant. 

Above the homosphere lies the heterosphere (100-500 km), where the primary 
transport mechanism is molecular diffusion. Very little vertical mixing occurs in this region, and 
the molecules tend to separate based on differences in molecular weight. 

Finally, the exosphere (above 500 km) identifies the layer where collisions 
between atoms become extremely infrequent. Because of the infrequency of collisions, particles 
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have high kinetic energies, and the escape of molecules from the atmosphere becomes 
significant. 

2. Ionosphere 

The ionosphere (50-1000 km) encompasses the upper part of the homosphere, 
the heterosphere and the lower region of the exosphere. The relative lack of vertical mixing in 
this layer, combined with an increasing flux of solar photons results in the formation of a 
partially ionized plasma. Although electrically neutral as a whole, the atmosphere in this region 
is comprised largely of long-lived positive ions and free electrons. Due to the large difference in 
mass between an ion and electron (the ions can be 50,000 times as massive as the electrons), the 
electrons are much more mobile. Hence, it is primarily the free moving electrons which interact 
with electromagnetic waves in the ionosphere. By accurately measuring electron density profiles 
the electromagnetic properties of the ionosphere can be inferred and exploited tactically by the 
military. 

The primary source of free electrons in the atmosphere is photoionization, where 
an incident solar photon strikes a neutral atom or molecule, forming a positive ion and a free 
electron. The process can be described by the equation 

M + ho -» M* + e~ , (1) 

where M is any molecule (or atom) and hv is the photon energy. The rate of this reaction 
depends on the product of the concentration [M] and the flux of photons [hv]. At high altitudes 
there is a large flux of photons but very few particles, and at low altitudes the reverse is true, so 
that the electron production function peaks at some intermediate level. In 1931, Chapman 
devised a simplified theoretical model of this process by assuming the atmosphere was 
comprised of only one constituent. Despite this rather drastic simplification, the theory yields 



remarkably good results. Since each dominant atmospheric species has its own ionization 
potential and reaction rate, the actual atmosphere can be represented by the superposition of a 
number of Chapman layers. A typical Chapman profile is shown in Figure 2. 



particular wavelength varies with altitude. In the EUV high energy, short wavelength photons 
generally penetrate farther into the atmosphere than long wavelength photons. The relative 
densities of different molecules also vary with altitude, and these effects combine to divide the 
ionosphere into four layers (D, E, FI and F2) based upon relative maxima in the electron density 
profile. These layers vary with time of day, geographic location and solar activity. 






The D layer (50-90 km) is the lowest region of the ionosphere and is produced 
by the penetrating component of the incident radiation, primarily the ionization of nitric oxide by 
the solar Lyman-a emission at a wavelength of 1216-A. At night, when there is no incident 
radiation, the electrons and ions quickly recombine and the layer disappears, except at auroral 
and polar latitudes. 

The next region, the E layer (90-140 km), is formed predominantly from the 
ionization of molecular oxygen by incident ultraviolet light and soft X-rays from the sun. Like 
the D layer, the E layer disappears at night except at latitudes above 60° where particle 
bombardment tends to maintain the ionization. 

The final region is broken down into two layers, the FI (140-200 km) and the F2 
(above 200 km). The FI layer is really a ledge in the electron density profile at the bottom of the 
F2 layer. It occurs only during daylight hours and is most pronounced during the summer months 
and at times of high solar activity (sunspot maximum). The presence of the FI ledge can be 
attributed to the ionization of molecular oxygen. The F2 layer has the highest electron density 
and is present during the day and at night. As with the FI layer, the F2 predominates during the 
daylight hours and summer months. However, unlike the FI region, the F2 persists at night 
because the density of particles is so low that collisions leading to the recombination of electrons 
and ions are extremely rare. A plot of typical electron densities versus altitude for both daytime 
and nighttime conditions is shown in Figure 3. 
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3. Electromagnetic Wave Propagation 

The ionospheric influence on radio wave propagation is significant and is the 
basis for a great deal of study. The radio frequency (RF) atmospheric window ranges from 10 
MHz-20 GHz. Due to their interaction with electrons, radio wave frequencies must be higher 
than the “plasma frequency” in order to pass through the atmosphere and reach a satellite or 
probe. The plasma frequency is the natural oscillation that electrons have in a magnetic field. For 
OTHR and “skip” communications, the frequency must be below the plasma frequency so that 
reflections of radio waves will occur. An ability to distinguish between propagation modes under 
a variety of circumstances establishes the requirement for an accurate electron density “map” of 
the ionosphere. 
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III. THEORY 


This section introduces the basic theory required to understand the method, results, 
discussion and conclusions regarding the AMI. First a general overview of interference and 
diffraction theory is presented. For a more detailed presentation on interference and diffraction 
theory refer to any standard text on basic physics, electromagnetic theory or optics, such as 
Hecht (1987) or Mailer (1988). This is followed by a description of the AMI, solar line profiles, 
specifically the oxygen 1304-A emission, solar occultation and a short description of the 
mathematical equations used to model the line profile of interest. 

A. NATURE OF LIGHT 

The work of Maxwell and others has shown that light is electromagnetic in nature, and 
classical electrodynamics points to the idea of a continuous transfer of energy by way of 
electromagnetic waves. Quantum electrodynamics, on the other hand, describes electromagnetic 
interactions and the transport of energy in terms of massless elementary particles known as 
photons, which are localized, propagating quanta of electromagnetic energy. This dual nature of 
light is evidenced by the fact that it propagates through space in a wave-like fashion and yet can 
display particle-like behavior during emission and absorption processes. 

If the wavelength of light is small in comparison to the size of an instrument’s optics, the 
techniques of geometrical optics may be used. A more precise treatment, applicable when the 
dimensions of the apparatus are small, is that of physical optics. In physical optics the dominant 
property of light is its wave nature. 



B. INTERFERENCE 


Interference is a consequence of the wave nature of light and can be described as the 
phenomena that occurs when light rays combine. The intensity of the resultant wave may be zero 
at certain positions, creating dark fringes and an interference pattern. 

Light waves can be considered as vector quantities comprised of two transverse vector 
components propagating together, the electric field and the magnetic field. The irradiance I, 
which is proportional to the time average of the magnitude of E squared (£ • 2T = (E 2 )), is of 
primary interest for most optics applications because it, or some quantity generated by it, is 
typically sensed by the optical detectors. 

Consider two coherent plane waves of equal magnitude, emanating from the same source 
and represented by their vector components Ei and E 2 . The complex representation of these 
waves traveling in the positive x direction are 

E t = Ae-’ 1 ^ 0 (2) 

E 2 =Ae- Hb ^ ,H) , (3) 

where A is the amplitude, co is the wave frequency, k is the propagation vector (k=2jtA.) and <(> is 

the phase shift between the two waves. A plane of constant phase that is perpendicular to the 
direction of propagation is called a wavefront. According to the principle of superposition, the 
electric field intensity E, at a point in space arising from the separate fields Ej and E 2i is given by 
E = E y +E 2 = Ae' Kkx ^‘\l + e~‘* ). (4) 

The last term may be rewritten as 
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1 + e~" 


(5) 



and the irradiance becomes 


/=£.£• =4^co S ^|), (6) 

where the term 2zl cosj^—J is called the amplitude factor. It is evident from equation (6) that the 

time dependence is automatically eliminated when calculating the irradiance, and I depends only 
on the square of the amplitude factor. 

In many experiments, relative intensity measurements are recorded, and the irradiance 
can be described by 



where I 0 is the reference intensity. Constructive interference occurs when the absolute value of 
the amplitude factor is a maximum, and destructive interference takes place when this factor is 
zero. In other words, constructive interference takes place where §=2mn and destructive 
interference where <j>=(2m+l)jt, for any integer m. 

Interference of light waves from different sources is the periodic enhancement and 
reduction of the irradiance of the superimposed waves. It is caused either by the propagation path 
length difference or an explicit phase shift between the waves or both. Interference theory can be 
extended to any number of superimposed waves, not just the two wave interference example 
used in the above derivation. 

The interference of light waves can be produced by splitting a single wavefront into two 
waves which travel different paths. Double-slit interference is an example of this and was first 
demonstrated by Young in 1802. Monochromatic, collimated light is allowed to strike an opaque 
screen which contains two small slits a short distance apart, and the resultant pattern is viewed 
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on a second screen some distance away from the first screen. An illustration of this is depicted in 
Figure 4. In this figure, y is the distance between the slits, x the distance from the slits to the 
detector, 0 the angle between the centerline and the observation point P and d the difference in 
distance between optical paths. 



Figure 4 Double-Slit Interference 


Since the light is monochromatic, the propagation vector for the wave leaving each slit 
has the same magnitude, 2%IX. The electric field at P is the superposition of the fields coming 
from each slit at that point, and the irradiance can be described by Equation (7). Assuming the 
detector is a large distance from the slit, the path difference can be approximated by dnysinG, 


and the phase shift between the separate waves becomes 










( 8 ) 



Substitution of Equation (8) into Equation (7) gives an irradiance at P of 

,,; iCOS .(SSl). (9) 

From Equation (9) it is evident that I is a maximum when the optical path difference is 
an integer multiple of X and a minimum when the path difference is a half integer multiple of X. 
By restricting the observation point on the detector to a small area around the optical axis such 
that the off-axis distance z«x, sin0 can be approximated by 0 (sin0«0«z/x). Equation (9) then 
reduces to 

/ = / 0 cos 2 ^^j, (10) 

and the detection screen will display equally spaced bright fringes separated by a distance of 

The intensity of the double-slit interference pattern will decrease as z increases for two 
reasons. The first deals with the initial assumption of infinitesimally narrow slits radiating 
uniformly in all directions. Due to diffraction effects, however, light does not radiate uniformly 
but decreases to zero for large values of 0. Diffraction is discussed in more detail in the next 
section. 

The second reason the intensity of the pattern falls off to zero is the concept of coherence 
length. Light does not travel through a medium as a perfect, continuous sinusoid but instead as 
wave groups that appear approximately as sinusoids for some specific length, defined as the 
coherence length. This is described as a photon wavetrain. The importance of coherence lies in 
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the fact that for interference to take place, paired groups of similar wavetrains must arrive almost 
simultaneously at P, If the path length difference is greater than the coherence length, wavetrains 
do not arrive within the coherence time, and interference does not occur. For light sources such 
as lasers the coherence length can be long (~1 m), and the interference pattern is easy to produce. 
For most natural light sources such as incandescent lamps and discharge tubes, which show no 
particular coordination between emitting atoms, the coherence length is short, and the 
interference pattern can only be generated with very small optical path differences. That is the 
reason interference patterns are much easier to produce with laboratory lasers. 

C. DIFFRACTION 

Diffraction is a general characteristic of wave phenomena and occurs whenever a portion 
of a wavefront (sound, matter or light) is obstructed in some way. The distinction between 
interference and diffraction is somewhat artificial, but it is customary to speak of interference 
when considering the superposition of a few waves, and diffraction when referring to a large 
number of waves. The mathematical analysis of this phenomena is generally treated by one of 
two approximations: Fresnel, or near-field diffraction, and Fraunhofer, or far-field diffraction. 
Fresnel diffraction occurs when the plane of observation is in close proximity to the obstruction, 
and Fraunhofer diffraction when the observation plane is a large distance from the obstruction. 
This thesis deals with far-field diffraction of light waves, so only Fraunhofer diffraction of light 
waves will be discussed. 

Diffraction effects are a consequence of Huygens-Fresnel principle, which states that 
every unobstructed point of a wavefront serves as a source of spherical secondary wavelets. The 
amplitude of the field at any point beyond the obstruction is the superposition of these wavelets, 
in both amplitude and phase. An illustration of Huygens-Fresnel principle is depicted in Figure 
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5. As shown in the diagram, light traveling through the central portion of the aperture continues 
to propagate in the same direction, but light passing by the edges of the aperture bends around 
the obstruction and propagates in a different direction. Near the aperture, the resulting waves 
appear as plane waves with slightly curved edges. When viewed farther away, however, the 
wavefronts become spherical. 

If the single aperture of Figure 5 is replaced by a number of slits, the interference 
phenomena described in the previous section are superimposed on the diffraction pattern from 
each of the slits. The resulting pattern is a series of alternately spaced light and dark fringes. This 
is pictured in Figure 6. 
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The intensity of the bright fringes is a maximum at the center of the pattern and drops off 
to zero as the distance from the optical axis increases. The overall effect is to create a number of 
different directions where constructive interference occurs. These are termed diffraction orders 
and are designated by m—1, 0, 1, etc., depending on their direction. 

The equation describing the relationships of Figure 6 is called the grating equation and is 
defined as 

/wA. = asin0 m , (12) 

where m is the diffraction order, X is the wavelength, a is the distance between slits and 9 m the 
deflection angle measured from the grating normal. For a source with a broad continuous 




spectrum, the m=0 order image corresponds to the unaltered, 0 o =O, white-light view of the 
source. Since Equation (12) is dependent on X , any value of m*0 results in a number of colored 
beams diffracted at slightly different angles. 

D. THE DIFFRACTION GRATING 

A diffraction grating is a repetitive array of diffracting elements, such as slits, that alters 
the phase or amplitude of a wave. The effect of a diffraction grating is very similar to a multiple 
slit arrangement, and the properties of the grating can be described using a modified version of 
the grating equation. Figure 7 shows a diagram of a typical reflecting diffraction grating. 



Figure 7 


A Diffraction Grating 




For this arrangement, the grating equation becomes 

mk = o(sin0 m -sinG,), (13) 

where a and m are defined as in Equation (12), 0j is the angle of incidence and 0 m is the 
diffraction angle for that order. 0 m is measured clockwise from the grating normal and 0j is 
measured counterclockwise from the normal. 

E. THE ALL-REFLECTION M3CHELSON INTERFEROMETER 

The All-Reflection Michelson Interferometer (AMI) is a modified version of a 
conventional Michelson interferometer which utilizes a blazed diffraction grating vice a partially 
silvered mirror to split and to redirect the incident light rays. A schematic diagram of a 
conventional Michelson interferometer is shown in Figure 8. 

As shown in the diagram, collimated light strikes a beam splitter positioned at angle of 
45° to the incident radiation. The beam splitter separates the light into two rays, forcing the 
resultant plane waves to travel in different directions. One beam is reflected from a fixed mirror 
back to the splitter while the other beam is reflected off a movable mirror. The two rays are then 
recombined by the beam splitter and directed towards the detector. Adjusting the movable mirror 
will change the distance one of the beams travels, changing the interference pattern at the 
detector. If the path length difference between the two waves is zero or a multiple of the 
wavelength, constructive interference will occur. If the path length difference is a multiple of one 
half a wavelength, destructive interference will take place. Translating the movable mirror in one 
direction will cause the intensity at the detector to alternate from bright to dark. 

It was not possible to use a conventional interferometer design for our objectives 
because the emissions of interest fall into the EUV and FUV portion of the electromagnetic 
spectrum, and most optical elements are opaque to light in this region. Several all-reflection 
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interferometers utilizing diffraction gratings as beam splitting devices have been built and 
evaluated (Roesler and Harlander, 1990; Fonck et al., 1978; Kruger et al., 1973). 



The prototype instrument examined in this thesis was first developed by Wallace (1992). 
The AMI is shown in Figure 9 and consists of a beam collimator, a plane diffraction grating, two 
plane mirrors and a detector. The collimator can contain a number of elements but is usually 
comprised of an aperture (pinhole or slit) and an off-axis parabolic mirror. The aperture is placed 
at the focus of the parabolic mirror so that light reflected from the mirror is collimated. In the 
diagram, the minus one order beam is represented by solid arrows and the zero order beam is 
shown as dotted arrows. The ray paths are also numbered to indicate the sequence of reflections 
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from when beams first hit the grating to when they recombine at the detector. Light striking the 
blazed grating is split along two separate paths corresponding to the zero and minus one diffraction 
orders. Plane mirrors are positioned to reflect the incident beams directly back to the diffraction 
grating. As a matter of convenience, the plane mirrors are designated by the diffraction order of the 
light they reflect back to the grating; the minus one order mirror reflects the minus one order beam 
and the zero order mirror the zero order beam. To ensure that the optical path lengths for both rays 
are equal, the zero order mirror is mounted on an adjustable base. The beams are again diffracted 
by the grating with the exception that the diffraction orders are reversed. In this way, the light 
recombines along a common path from the grating to the detector and creates an interference 
pattern. 



22 




Recall from Equation (13) that the direction in which a wavefront propagates after 
striking the grating depends on the radiation wavelength ( X ), the diffraction order (m) and the 

angle of incidence (0,). The dependence on wavelength, in particular, allows for unique 
characteristics and forms the basis of operation of the AMI. Figures 10, 11 and 12 illustrate the 
effects of wavelength and diffraction order on the plane waves, and the general operating 
principles of the AMI. 

Figure 10 details the initial reflection of the plane waves from the diffraction grating. As 
in Figure (9), the minus one order plane waves are shown as solid lines while the zero order 
waves are pictured as dashed lines. 
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The zero order beam shows no angular dispersion because the angle of incidence equals 
the angle of reflection when m=0. The minus one order beam, on the other hand, experiences an 
angular shift of A0 due to its diffraction order (m— 1) and X. 

Figure 11 shows the effect the plane mirrors have on the propagating wavefronts. The 
zero order beam approaches perpendicular to the mirror which means it is reflected 
perpendicular to the mirror. The minus one order beam travels toward the mirror with an angular 
dispersion of A0; therefore, it reflects off the mirror with the same A0. The mirror shifts the A0 
to the opposite side of the normal, but it does not add any additional angular shift. 
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Figure 12 illustrates the second reflection of the beams from the grating and the 
recombination of the wavefronts at the detector. For the second reflection the diffraction orders 
are reversed. The minus one order beam undergoes an m=0 reflection while the zero order beam 
experiences an m=-l reflection. This means that the zero order beam, not the minus one order 
beam, exhibits an angular shift of Aft after reflecting from the grating. 



As Figures 10, 11 and 12 show, each ray undergoes an angular shift of Aft at some point 
along its optical path length. Each ray is shifted by the same amount, but in opposite directions, 
so that the total angular dispersion at the detector is 2A0 (where 0 depends on X). This is 
illustrated in Figure 12 and the enlarged view at the detector of Figure 13. 
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By adjusting the mirrors to exactly offset the angular shift of the wavefronts for 


particular wavelength, Ao, the rays will arrive at the detector parallel to each other and interfere 
constructively at all points. In other words, the light at the detector will appear identical to the 
incident light. When the AMI is adjusted for a particular wavelength in this manner, the 
wavelength is referred to as the “tuned” wavelength (Kq). 



Figure (14) represents a specialized case of Figure (13) where the mirrors of the AMI 
have been adjusted to produce one fringe cycle across the detector. The lowest spatial frequency 
occurs when one bright fringe and one dark fringe strike opposite ends of the detector with 
length L. The diagram of Figure (14) and Equation (13) can be used to calculate the resolution 
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Once 0 has been determined for the tuned condition, Equation (13) can be manipulated 
to determine the relationship between changes in wavelength and changes in angular dispersion. 
For small incremental changes, 80 can be approximated by A0, and the derivative of Equation 
(13) becomes 

AX = a cos(0 )A0 . (16) 

Performing a trigonometric identity substitution for the cosine term and replacing the resulting 
sine term with Equation (13) yields 

AX = aAQ Vl — sin 2 (9 ) = aAQ , (17) 

which is valid for any value of A0. 

Using the small angle approximation (tan(A0)«A0), A0 from Figure (14) can be 
expressed as 

A0=M1 (18) 

With this relationship between A0, Xq and L, Equation (17) becomes 



where A)^,, n is the limit of resolution for the instrument. Using this method, Wallace (1992) 
calculated a typical resolution of the AMI to be approximately 4 mA. 

For the AMI, which uses a blazed diffraction grating. Equation (13) can also be written 
as 


pmX = sin0 m - sinG,, 


(20) 



with p representing the ruling density of the grating in lines/mm (p=l/a). Cleary et al. (1992) 
have shown that judicious selection of p with respect to X will maximize the light incident on the 
detector by eliminating all reflections of unwanted diffraction orders. Using the conditions that 
9i and all higher order diffraction orders are greater than 90° and 0. 2 and all lower order 
diffraction orders are less than or equal -90°, Equation (20) can be solved for p to give 
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( 21 ) 


where p min is the minimum acceptable ruling density for the AMI. Using Equation (21), Wallace 
(1992) has shown that p m[n =6000 lines/mm is desirable for the oxygen 1304-A emission lines. 

F. THE SOLAR OXYGEN 1304-A TRIPLET LINE PROFILES 

The solar spectrum is dominated by strong emission lines below approximately 1500-A. 
Although the strongest of these lines are attributed to atomic hydrogen and helium, there is also a 
strong atomic oxygen emission at 1304-A. The brightness of this emission is due to the 
accidental resonance between the Lyman-p line of hydrogen at 1026-A and a ground state 
transition of atomic oxygen. This resonance is known as the “Bowen mechanism” and it 
produces three oxygen emission lines at 1302.17, 1304.86 and 1306.02-A. Since the solar 
atmosphere is extremely opaque at these wavelengths, all three lines display a pronounced self 
reversal at the center of the line profiles. The 1304.86-A emission line was examined in this 
thesis, but the other two emission lines are similar and can be modeled in the same fashion. 

The Doppler-broadened 1304.86-A line profile can be approximated as the sum of two 
equal and offset Gaussian curves according to the equation 


2Vtt x d 


( 22 ) 
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where rcF(x) is the flux in ergs cm' 2 sec' 1 , itF is the line-integrated flux in ergs cm' 2 sec' 1 , is a 
measure of the line dispersion in km/sec and x ofT is a measure of the offset of the emission peaks 
from the center of the line in km/sec. This method was developed by Gladstone (1992) who 
generated an empirical fit to data obtained from the Solar Maximum Mission (SMM). Equation 
(22) is expressed in terms of flux and velocities. For the purposes of this thesis, it is more 
convenient to speak in terms of relative intensity and wavelengths. Expressed in these units, 
Equation (22) becomes 



where I 0 is the maximum intensity and the offset and dispersion coefficients are as defined above 
(in units of A). A model of the 1304.86-A emission line using the values X of f=.0478-A and 
A dl5 =.0478-A is shown in Figure 15. The Interactive Data Language (IDL) procedures used to 
generate the profile are listed in Appendix A. 
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Figure 15 Solar Oxygen 1304.86-A Emission Line Profile 


G. THE TERRESTRIAL OXYGEN 1304-A EMISSION 

The oxygen 1304-A triplet emission is also a prominent feature of the terrestrial airglow. 
The terrestrial emission is caused by solar resonant scattering and photoelectron impact 
excitation of atomic oxygen in the thermosphere. The terrestrial emission differs from the solar 
emission in that it displays a narrower lineshape and the intensity is much weaker (the solar 
emission is over 10 4 times as bright as the terrestrial emission). 

The shape of the terrestrial oxygen 1304-A line profile is dominated by two line 


broadening mechanisms, Doppler broadening and natural broadening. Of the two, Doppler 




broadening is probably the easiest to understand. Atoms in a gas are always in random thermal 
motion, and the average speed of the atoms is dependent on the amount of thermal kinetic energy 
they possess. From the point of view of a stationary observer, this means that at any instant the 
atoms will be traveling at different velocities, some towards the observer and some away from 
the observer. As a result, the frequency of emission or absorption of each atom in its own frame 
corresponds to a different frequency from that for an observer. Each atom has its own Doppler 
shift which spreads the emission or absorption line out (although the total strength remains 
constant). The broadening of the line can be described by a simple Gaussian line profile. The 
relationship between temperature and the Doppler width of the spectral line is given by 


^Do PP ,er = 


Vj_ f2kT 

d M ' 


(24) 


where Av is the thermal Doppler width (in units of frequency), v 0 is the line center frequency, c 


is the speed of light, k is Boltzmann’s constant, T is the temperature and M is the atomic mass. 
In terms of Doppler wavelengths. Equation (24) becomes 


A , _h. He 

A Vpto- c V M ' 


(25) 


The line-center cross section for each al 


= -fn- 
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(26) 
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for the case of Doppler broadening. Replacing the constants of Equation (26) with numerical 
values gives 


a Vo = 1.6*10-% JI77/ 12 , 


(27) 
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with Xq in A and T in K. A is the atomic weight of the atom and f| 2 is the oscillator strength. 
Oscillator strength is defined as the ratio of the radiative decay rate of a real atomic transition to 

the decay rate of a classical electron oscillator at the same frequency. Since it is a ratio of rates, 
the oscillator strength is dimensionless. 

Natural broadening is a result of the Heisenberg uncertainty principle which relates the 
energy and lifetime of a system. In equation form the uncertainty principle is 

AEAt >h/2, (28) 

where AE is the uncertainty in the energy of the system. At is the time interval characteristic of 
the rate of change of the system and ti=h/2n. If the lifetime of an atomic state is long (approx. 
10~ 3 sec) than the line profile will be relatively narrow, whereas atomic states that are short-lived 
(approx. 10 s sec) display broad lineshapes. For natural broadening, the widening of the spectral 
line can be described by a Lorentzian line profde. 

Quite often, as in the case of the terrestrial 1304-A emission, atoms exhibit both natural 
and Doppler line broadening effects. In these situations, the lineshape can be described by the 
Voigt function, which is a convolution of the Lorentz profile over the Doppler profile. The 
terrestrial 1304.86-A emission line profile can then be expressed as 

I Terrestrial = hi*™ dY'™ ^ H(a,U), (29) 

where H(a,u) is the Voigt function, and the Voigt parameters are defined as 


r 

4?rAv c ’ 

(30) 
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U Av z) 

(31) 

r=y„ +y, +2v c . 0/ . 

(32) 
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The parameters y u and yi are the widths of the upper and lower states involved in the transition 
and v col is the collision frequency. 

A more detailed explanation of the Voigt function, and line broadening mechanisms in 
general, can be found in Rybicki and Lightman (1979). Figure (16) shows a comparison of 
Doppler and Voigt profiles for the terrestrial oxygen 1304.86-A emission line. The IDL 
procedures used to model the profiles are listed in Appendix A. The terrestrial lines were 
generated using a temperature of 300 K and a Voigt parameter of a=1.108xl0* 2 . As the diagram 
shows, the shape of the profiles are very similar near the center, but the Voigt curve widens out 
significantly near the edges while the Doppler profile remains fairly narrow. This indicates that 
Doppler broadening dominates near the center of the line profile while natural broadening and 
the Lorentzian lineshape is predominate in the regions away from line center. 
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Figure 16 Comparison of Doppler and Voigt Line Profiles 


H. SOLAR OCCULTATION 

Observations of the solar oxygen 1304.86-A emission line from within the atmosphere 
will be attenuated to some degree by the intervening column of oxygen, and the amount of 
absorption is directly proportional to the column density. The larger the column of oxygen 
between the observer and the solar emission, the greater the attenuation of the solar line. The 
equation describing the resulting lineshape is given by 

/ 0iKWd (A.) = (33) 
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where N is the column density (cm' 2 ) and cr (A.), the resonant scattering cross-section (cm 2 ), is 
given by 

cr(A) =cr 0 Av o II(a,u). (34) 

Since a is also dependent on temperature, there are really two geophysical variables m 
Equation (33), N and T. Figures (17) and (18) show calculations of the variation of Iobserved with 
changing column densities. These calculations assumed an isothermal atmospheric temperature 
of 300 K and a Voigt parameter a of 1.108xl0' 2 . Figure (17) plots four curves, the unattenuated 
solar line profile and three other curves corresponding to column densities of lxlO 17 , lxlO 18 and 
lxlO 19 cm' 2 . Since these column densities are representative of the actual column densities in the 
atmosphere, Figure (17) is a fairly accurate simulation of observed oxygen 1304.86-A emission 
line profiles. 
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Figure 17 Simulated Oxygen 1304.86-A Line Profile versus Column Density 
(N=lxl0 17 , 1 xlO 18 and lxlO 19 ) 

Figure (18) plots the simulated line profile for column densities of lxlO 21 , lxlO 22 and 
lxlO 23 cm' 2 and shows that at high column densities, the solar 1304.86-A is almost completely 
attenuated by the atmosphere. These figures also indicate that any measurement technique of the 
observed oxygen 1304.86-A line profile will be sensitive to the column of oxygen over the entire 
altitude range of the thermosphere. The IDL procedure used to model the observed intensity is 
shown in Appendix A. 
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Simulated Oxygen 1304.86-A Line Profile versus Column Density 
(N=lx10 21 , lxlO 22 and lxlO 23 ) 


The profiles of Figure (19) illustrate the expected temperature dependence of an 
occultation measurement. The solid lines of Figure (19) depict the column densities of Figure 
(17) for a temperature of 300 K while the dashed lines plot the same column densities for a 

temperature of 1000 K. 
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Figure 19 


Simulated 1304.86-A Line Profile versus Temperature 


Figure (20) represents an enlarged view of Figure (19) around the center wavelength. 
Figures (19) and (20) show that higher temperatures broaden the absorption of the simulated line 
profile around the center at the same time as decreasing the depth of the attenuation. Although 
the total amount of attenuation remains fairly constant, the shape of the line profile is definitely 
affected by the temperature and must be accounted for when analyzing the iineshapes. 
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Figure 20 Enlarged View of Line Profile versus Temperature 

Figure (21) is an illustration of the curve of growth for the total simulated oxygen 
S-A emission at a variety of temperatures. The IDL procedure used to plot the growth 


in Appendix A. 





Simulations of AMI performance by Carlson (1992) using the Dart ray tracing program 
show that AMI is an excellent instrument for observing the oxygen 1304.86-A emission line. 
Once the line profile and temperature at various altitudes have been measured, the amount of 
atomic oxygen in the thermosphere can be determined using methods such as that developed by 
Meier (1992), which utilizes Discrete Inverse Theory to retrieve the oxygen concentration from 
solar occultation data. 
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IV. EXPERIMENT 


This section describes the experimental procedures used to setup and test the prototype 
All-Reflection Michelson Interferometer. 

A. SETUP OF AMI 

The initial setup and testing of a prototype AMI was done by Risley (1992). He built an 
instrument similar to that described in Figure (9), with the exception that he inserted a lens 
between the light source and pinhole aperture in order to focus as much light as possible on the 
pinhole opening. After construction was complete, he proceeded to test the instrument using a 
number of different light sources. Unfortunately, he achieved only limited success due, 
primarily, to an inadequate detector. For a more detailed treatment of the original construction 
and prototype testing see Risley (1992). 

The first step in the continuation of prototype testing was to increase the intensity of 
light reaching the detector. This was done by comparing the original grating (p=1200 lines/mm) 
to a sinusoidal diffraction grating with a ruling density of 2400 lines/mm. It was hoped that a 
sinusoidal grating with a higher ruling density would increase the intensity of light in the zero 
and minus one order beams. Measurements of beam intensity were made with both gratings 
using two different light sources, a He-Ne laser and a mercury vapor lamp. Figure 22 illustrates 
how the beams were designated and where readings were taken. Tables 1 and 2 are a compilation 
of the results. Table 1 shows the measurements in terms of light intensity (watts) while Table 2 
depicts the results in terms of percentages. The m=(0,-l) and m=(-l,0) percentages of Table 2 are 
based on the intensities from the m=0 and m— 1 beam, not the initial beam. For example, the 
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TABLE 2 


PERCENTAGE COMPARISON OF GRATINGS 
(Standard with p=1200. Sinusoid with p=2400) 


Reflection 

He-Ne Laser (%) 

Mercury Lamp (%) j 


1200 

2400 

1200 

2400 

Initial 





m=0 

15.2 

69.5 

21.0 

70.8 

m=-l 

66.7 

16.8 

35.4 

9.2 

m=(0,-l) 

19.3 

13.1 “1 

9.4 

6.4 

m=(-l,0) 

5.7 

57.1 

6.9 

53.4 

Combined 

6.7 

19.0 

4.4 

9.6 


The results show that the sinusoidal grating achieved better results for both light sources 
than the standard grating. In fact, the sinusoidal grating showed a twofold increase in the light 
reaching the detector. The light intensity for the He-Ne laser went from 6.7% to 19.0% while the 
intensity for the mercury lamp went from 4.4% to 9.6%. Interestingly enough, the sinusoid 
grating distributed the light differently between the two diffraction orders. The sinusoidal grating 
sent more of the light to the zero order reflection while the standard grating sent the majority of 
light to the minus one order. 

The next improvement made was to mount the optical workbench on a vibration 
isolation base. The Micro-g Series 62 post isolation system consisted of four supports bolted 
together with tiebars to comprise a structurally sound rectangular base. Each support contained 
an air piston and leveling valve for final adjustments. A pressurized nitrogen gas tank provided 
air to the pistons and main pressure was maintained at approximately 60 psi with a regulating 

The third refinement, shown in Figure 9, was to mount the zero order mirror on an 
adjustable base to insure the path lengths for both rays were equal (within the coherence length 
of the light source). 
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The final, and perhaps most important, modification to the original prototype was to 
replace the Panasonic CCD video camera with a solid state Electrim EDC-1000HR integrating 
CCD camera. The EDC-1000HR consisted of a camera head assembly , a computer interface 
card and connecting cable and was controlled directly from the computer keyboard. The sensing 
area of the detector was 8.67 mm wide by 6.59 mm high and the pixel array (resolution) was 753 
horizontally by 244 vertically. One of the main advantages of the EDC-1000HR camera was that 
images on the computer screen could be saved directly to hard drive or floppy disk as TIFF files 
for later manipulation and analysis. Other important features included the ability to control 
exposure time from 1 msec to 1 minute and to subtract background light from the picture. 

B. TESTING WITH MERCURY VAPOR LAMP 

The first step in testing the interferometer design was to create an interference pattern 
using a light source with a short coherence length. A low pressure mercury discharge tube 
mounted in an Oriel model 65160 lamp, such as that used by Risley (1992), was chosen. A band¬ 
pass filter was placed across the opening of the lamp to filter out all emission lines from the 
discharge tube except the green line of interest at a wavelength of 5461 -A. 

The setup was exactly the same as Risley’s (1992). A convex lens was placed between 
the lamp and a 250 pm pinhole and the pinhole was placed at the focus of the lens to maximize 
the light throughput. The Oriel Model 45347 off-axis parabolic mirror was focused as before, 
using a small telescope focused at infinity. Since the sinusoidal diffraction grating increased the 
light intensity of the beams, the zero and minus one order mirrors were aligned by observing the 
reflections from the grating and placing the mirrors visually in the beams so that the light 
bounced off the mirrors and returned to the same spot on the diffraction grating. 
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Once the initial setup and alignment was complete, the computer camera was turned on 
and actual testing commenced. It took only a few minutes to adjust the mirrors and establish an 
interference pattern. Figure 23 is a picture of the first interference pattern achieved with the 
mercury vapor lamp at 5461-A. 


Figure 23 Initial Interference Pattern for Mercury Vapor Lamp 
The exposure time for the interference pattern was 1 sec and the image was recorded 
using the HWINCAM software provided with the computer camera. The HWINCAM program 
only provided a pixel resolution of 192 (H) x 180 (V) and was used only for the first picture. All 
subsequent images were recorded using the HVESA program (also supplied by Electrim) which 
provided the maximum resolution of 753 (H) x 244 (V). 
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Once the pattern was established, fine adjustments to the system were made to improve 
the contrast and center the pattern. The minus one order mirror had two adjustment screws, one 
to control the vertical angle and another the horizontal angle of the mirror. Adjusting the vertical 
screw affected the tilt of the interference lines. If the screw was turned far enough in one 
direction, the fringes rotated a full 360 degrees. The adjustment screw was extremely sensitive 
and it proved difficult to get the lines perfectly vertical. The horizontal screw controlled the 
spatial frequency of the interference pattern. Rotating the screw clockwise or counterclockwise 
moved the mirror angle closer to, or farther away from, the tuned wavelength and affected the 
number of interference fringes on the screen. As the wavelength approached Xq the number of 
fringes on the screen decreased. Moving the wavelength away from hj had the opposite effect 
and increasing the mirror angle such that the wavelength was too far from Xq caused the fringes 
to disappear altogether. 

Another important adjustment was the position of the zero order mirror from the grating. 
When the pattern was first established, the brightest fringes were near the right edge of the 
screen. Moving the mirror closer to the grating caused the bright fringes to move towards the 
center of the screen. When the distance from the zero order mirror to the diffraction grating 
equaled the distance from the minus one order mirror to the grating, the brightest fringes were in 
the center of the screen and the interference pattern displayed the greatest definition. 

Experimentation with screw adjustments and settings showed that the AMI was also very 
sensitive to vibrations and air currents around the mirrors. Even with the workbench resting on a 
cushion of air, footsteps or movement around the table caused the pattern on the screen to 
oscillate. The only way to create a stable pattern was to make adjustments and then stop all 
motion for a few minutes to give the interference pattern on the screen time to settle. 
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Since the AMI was working as predicted, all that remained was to record and analyze 
data for a number of sources. Three IDL procedures were written to analyze data recorded by the 
camera. The first procedure, named IMAGE, read information from an interference pattern 
(stored as a TIFF file from the EDC HVESA program) and displayed the image on the screen. 
The second procedure, called PROFILE, vertically summed the intensities of bright and dark 
fringes from the interference pattern. The net result was an interference profile (753x1 array) 
which showed the variation of light intensity as a function of fringe location on the detector. The 
final procedure, designated as SPECTRUM, took the Fourier transform of the interference 
profile to recover the initial spectrum. For the mercury emission at a wavelength of 5461-A, 
which approximates a monochromatic light source, the resulting spectrum was a single spike 
some distance from X<). A number of runs were conducted at different distances from the tuned 
wavelength and the following figures are indicative of the data collected. Figure 24 shows an 
interference pattern. Figure 25 an interference profile and Figure 26 the Fourier transform of the 
interference profile. 

The image of Figure 24 is very similar to that of Figure 23. Close inspection of both 
figures reveals a background pattern that resembles a fingerprint. Every possible step was taken 
in an effort to localize the source of the background pattern. First, each light beam was blocked 
separately to see if one of the optical elements was responsible for the pattern. It made no 
difference which beam was blocked, as long as any light reached the detector the background 
pattern was visible. All elements common to both beams were also inspected to no avail. Finally, 
the detector was disassembled and the glass covering the solid state surface was cleaned in an 
effort to remove the pattern. Unfortunately, the glass was permanently affixed to the detector 
surface which prevented the surface itself from being cleaned. Eliminating all other sources for 



the background pattern, it was surmised that the pattern was, indeed, an actual fingerprint which 
was deposited on the surface sor me during the manufacturing process. 



Figure 24 Interference Pattern for Mercury Vapor Lamp 
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The interference profile of Figure 25 is an excellent example of Fraunhofer diffraction. 
As expected, the number of peaks and valleys match the number of bright and dark fringes from 
Figure 24. Since the AMI measures wavelength relative to a particular tuned wavelength, the 
Fourier transform of the interference profile should produce a single peak at a position on the x- 
axis dependent on the number of bright fringes in the profile (number of fringes is a measure of 
spatial frequency). There are nine bright fringes in Figure 25 plus three or four not shown on the 
graph. Therefore, Figure 26 should display a single spike approximately 12 to 13 units from Xo 


SO 





(the zero point on the x-axis). Examination of Figure 26 reveals that it does indeed show a 
pronounced spike at that position. 



The large peak near the origin of Figure 26 represents the low frequency noise 
component of the signal. This peak was present in all of the spectrums recorded during the 
experiment. In order to distinguish the actual wavelength peak from the background noise, it was 
necessary to adjust the AMI so that the interference pattern consisted of at least four or five 
bright fringes. This had the effect of moving the peak away from the low frequency noise. 
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27, 28 and 29 


recorded using the S£ 


achromatic light source but at a 


higher spatial frequency than Figures 24, 25 and 26. The image of Figure 27 appears much 
cleaner and more defined than the image of Figure 24 because the background was subtracted 
from the interference pattern using the remainder feature of the EDC HVES A software. 



Figure 27 Interference Pattern for Mercury Vapor Lamp (High Spatial Frequency) 
Figure 28 shows 14 bright fringes, which means that, including the three or four fringes 
not shown, there should be a spike at approximately 17 or 18 units on the graph of Figure 29. 
Inspection of Figure 29 shows that the predicted spike does indeed appear. 
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Figure 28 Interference Profile for Mercury Vapor Lamp (High Spatial Frequency) 


Figures 26 and 29 show another interesting feature of the AMI. As spatial frequency 
increases and the distance from the tuned wavelength increases, the intensity of the peak 
decreases. 
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Figure 30 illustrates this effect. The three curves plotted on the graph were recorded 
using the mercury vapor lamp as the light source. The only difference between the curves is that 
they correspond to different spatial frequencies. Curve A was plotted for an interference pattern 
that contained 11 bright fringes. Curves B and C show interference patterns with 20 and 29 
bright fringes, respectively. Figure 30 clearly displays the relationship between spatial frequency 
and intensity. As frequency increases, the intensity of the peaks drop from 1600 for Curve A, to 
900 for Curve B, to 300 for Curve C. As a result, this becomes an important design consideration 
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when determining the tuned wavelength for the AMI. If Xq is too far from the oxygen 1304.86-A 
line, the intensity of the profi.e might be to low' to accurately reconstruct the spectrum. 



C. TESTING WITH SODIUM VAPOR LAMP 

The next phase in testing the interferometer design was to verify the performance of the 
instrument using a different, less intense light source, sodium was chosen because it exhibits a 
doublet emission at a wavelength of 5889-A. which is less intense than either the He-Ne laser 
tested by Risley (1992) or the mercury vapor lamp examined previously. The nature of the 
doublet emission (two distinct emission lines within 9-A of each other) also provided a unique 
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opportunity for testing the AMI. Not only could the performance of the interferometer be 
verified with low intensity sources, the resolving power of the instrument could also be 
examined. 

The resolution limit and bandwidth of the instrument for the sodium doublet emission 
can be determined from Equation (16) and Figure (14). Using 


(35) 


for the doublet emission, Equation (16) becomes 

(3i 

With L=8.7 mm, a=4.1667xl0' 4 mm/line (1/p), 0=28.5° and >.=5889-A for the sodium 
lamp, A>v rain =.2487-A. Since the doublet emission lines are separated by approximately 9-A, the 
resolving power is >7AA mm =72 pixel units. 

The bandwidth (BW) can be determined using the Nyquist sampling criteria, which 
requires a minimum of two samples per cycle. Therefore 


W = j^A>. mi „, (37) 

where N is the number of detector elements. With N=753, BW=93.6-A. 

After completing the numerical analysis, the mercury bulb in the Oriel lamp was 
replaced with a sodium vapor bulb. Since the sodium lamp operated at a slightly different 
wavelength, the instrument was realigned visually so that both beams again intersected at the 
detector. As before, once the mirrors were properly oriented it only took a few minutes to 
establish an interference pattern and fine tune the image. Figure 31 is a picture of the resulting 
interference pattern. 
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Figure 31 Interference Pattern for Sodium Vapor Lamp (Xq less than A., and X 2 ) 

The image of Figure 31 clearly shows two superimposed interference patterns. The 
pattern with four or five fairly wide bright fringes is the more dominant of the two and is a result 
of the emission line which is nearest Li- The other pattern, with approximately 55 narrow 
fringes, results from the emission line farthest away from /-o- For future reference, the two 
emission lines from the sodium vapor lamp are designated as A, t and X 2 . Since the mirrors were 
aligned visually and not to any particular wavelength, there was no way of telling from Figure 31 
where X! and X 2 lay in relation to Xq. A relatively simple method, and the one used throughout 
the remainder of this thesis, was to adjust the AMI so that one of the interference patterns 
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completely disappeared. When one of the patterns disappeared, the instrument was tuned for that 
exact wavelength. By trial and error, it was then a simple matter to determine if the tuned 
wavelength was outside or between the sodium doublet emission lines. The data of Figure 31 
was recorded for Xq less than X, and X 2 . 

The only difficulty which arose while testing the sodium lamp was vertically aligning 
both interference patterns simultaneously. When one pattern was aligned the other was slightly 
askew and vice versa. After a number of attempts, the effort to align both patterns at the same 
time was abandoned and the patterns were lined up as close to vertical as possible. The 
interference profile for the image of Figure 31 is shown in Figure 32. 



lOcfV' 200. 300 " 400 &Q0 jSOO' 

DETECTOR POSITION 


Figure 32 Interference Profile for Sodium Vapor Lamp (Xq less than X] and X 2 ) 
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As with Figure 31, the profile of Figure 32 also reflects the influence of two separate and 
disti ; wavelengths. The wavelength with the lower spatial frequency gives the curve its general 
shape while the wavelength with the higher spatial frequency gives the curve its jagged 
transitions. 


Figure 33 plots the Fourier transform of Figure 32. Contrary to our expectations, Figure 
33 does not show two emission peaks. 



For some unknown reason the second peak was not visible eventhough two distinct 
spatial frequencies can be discerned in Figures 31 and 32. Although the data from the mercury 
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lamp indicates that the intensity should decrease with increasing spatial frequency, at a distance 
of around 30 to 40 pixel units some type of peak should be visible. In order to determine the 
reason only one peak was apparent, a number of runs were conducted for different tuned 
wavelengths. The AMI was adjusted and data recorded for the following situations: Xo on one 
side of X, and X 2 , Xo between X., and X 2 and Xo on the other side of X., and X 2 . Figure 34 illustrates 
the case where Xo is approximately midway between X, and X 2 . 



Figure 34 Interference Pattern for Sodium Vapor Lamp (Xo between X t and X 2 ) 

As in the case of Figure 31, the pattern of Figure 34 also shows two superimposed 
interference patterns. The fringes of Figure 34, however, are approximately the same width since 
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Xq is almost the same distance from X| and k 2 . Figure 35 is a plot of the interference profile for 
the image of Figure 34. Again, two wavelengths are visible on the interference profile. 



The graph of Figure 36 displays the Fourier transform of Figure 35. This time there does 
seem to be a faint second peak around pixel 36. The peak is too faint however, to say positively 
that it is a product of the second emission line. In all, ten runs were conducted but the data for all 
ten runs was inconclusive. Since the instrument seemed to be working properly and the second 
emission line was visible in the interference patterns and interference profiles, there had to be 
some other reason it was not appearing in the Fourier transform. 
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One possible answer was the methodology used to analyze the raw data. The interference 
profiles were obtained by summing the intensities vertically (over all 244 vertical pixels) and 
then taking the Fast Fourier Transform (FFT). As noted previously though, the fringes were 
never perfectly vertical. When the fringes for k, were vertical, the fringes for were slightly 
tilted. To compensate for this, the array was summed by rows of two instead of 244 rows at once. 
In other words, the intensities for two rows were summed, then the FFT was taken for those two 
rows. After this step was completed, the resulting 122 FFT’s were summed to create the final 
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Fourier transform array. The IDL procedure used to sum the array over a certain number of rows 
at once is listed in Appendix A. 

Summing the array two rows at a time had a dramatic effect on the data. The second 
emission peak was clearly evident for all ten runs. The results of this analysis are listed in Table 
3. 


TABLE 3 SODIUM DOUBLET EMISSION ANALYSIS (Summed by Rows Of 2) 


Run 

*0 

Position 

X, Pixel 
Distance 

+ or- 

X 2 Pixel 
Distance 

Total 

Separation 

1 

A»|<X()<X.2 

14 

+ 

36 

50 

2 





50 

3 

Xo<^t 




51 

4 

A.| < 7.2<Xfl 

5 


55 

50 

5 

Xo<X.|<X<2 

4 


55 

51 

6 

X t <Xo<X 2 




50 

7 





50 

8 

Xi<Xo<X2 

18 

+ 

33 

51 

9 


25 

+ 

25 

50 


The results listed in Table 3 are fairly consistent. The data shows that, regardless of the 


location of the tuned wavelength, the instrument displays an average spacing between the 
doublet emission lines of 50 pixel units. In fact, the experimentally determined resolution of 50 
units compares favorably with the predicted resolving power of 72 units which was determined 
theoretically at the beginning of this section. This brings up an interesting point. When the 
resolving power was first determined theoretically, there was a mistake in Equation (36) and the 
resolving power was erroneously calculated as 36 units. After the runs were completed and the 
resolving power came out experimentally as 50 units, the equations were rederived, the algebra 
mistake was corrected, and the right theoretical value was determined to be 72 pixel units. 
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The findings of Table 3 bring up an interesting question. Namely, over how many rows 
should the data be summed to provide maximum definition between the peaks and the 
background noise? To answer this question the data from one particular run was examined. Run 
1 was chosen arbitrarily and the intensity profile was summed over 1 row, 4 rows, 16 rows, 64 
rows and 244 rows. The resulting five curves were then plotted and compared on the graph of 
Figure 37. Figure 37 shows that the definition between the second peak and the background 
noise improves steadily as the intensity distribution is summed up to about 64 rows at a time. At 
that point, the intensity of the second peak starts to decrease and it becomes harder to distinguish 
the peak from the background noise. Examination of other stored EDC files using this technique 
revealed the same general trend. Definition between the peak and background noise improved 
until the data was summed over approximately 32 to 64 rows at a time, then the definition started 
to decrease until the peak was indistinguishable from the background noise. The reason or 
reasons for this are not clearly understood at this time. 
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Figure 37 Comparison of Sodium Doublet Emission Profiles (Intensity Profile 
Summed Over 1,4, 16, 64 and 244 Rows) 

D. FINAL MODIFICATION OF AMI 

Validation of the prototype instrument as a tool to perform high resolution analysis of 
emission line profiles was finally complete. The only task that remained was to perform a slight 
modification to the prototype so that it reflected a design more suitable for operation in the 
vacuum ultraviolet. This was accomplished by removing the focusing lens and replacing the 
pinhole opening with a narrow slit. Removing the focusing lens between the light source and 
aperture enabled the instrument to collect any incident light within the geometric window of the 
slit-parabolic mirror configuration and gave the AMI a wider field of view. A slit was used in 
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place of the pinhole opening because it was easier to adjust the width of the slit and to control the 
amount of light striking the off-axis mirror. 

The AMI, with modifications complete, is pictured in Figures 38 and 39. Figure 38 
shows the front view of the instrument as seen from behind the detector. Figure 39 displays the 
rear view of the instrument viewed from behind the diffraction grating. The light source shown 
in Figures 38 and 39 is the Oriel lamp with mercury vapor bulb. The mercury lamp was used 
temporarily to adjust the various optical elements and see what effect, if any, the modifications 
had on AMI performance. 

In the front view of Figure 38, the zero order mirror (mounted on adjustable base) is the 
element in the left foreground of the picture. Directly behind the mirror is the diffraction grating 
mounted in its circular holder. To the right of the zero order mirror is the solid state detector with 
the computer cable attached to rear. Next to the detector is the minus one order mirror. The off- 
axis parabolic mirror with its polished surface can be seen in the background. The elements 
farthest to the right of the picture are the slit and, directly in front of it, the mercury lamp. 
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The slit itself can be seen in the rear view of Figure 39. The circular mirrors and detector 
opening are also visible. Another feature evident in the rear view of Figure 39 is the spacing 
between the diffraction grating and reflecting mirrors. The mirrors are situated as close to the 
grating as the base will allow. Tests showed that the distance between the grating and mirrors 
had no visible effect on AMI performance. In fact, the ability to fit the instrument in virtually 
any size container is one of its main advantages. 

Figure 40 shows an interference pattern recorded after final modifications to the AMI 
were completed. Figure 40 was taken with the mercury lamp and an exposure time of 500 msec. 
The interference pattern is not as distinct in Figure 40 as in previous figures because of the light 
intensity reaching the detector. Since the light source was directly in front of the slit, the amount 
of light reaching the detector was much higher and tended to wash out the picture. Even at half 
the exposure time (500 msec versus 1000 msec for other figures) it was difficult to get a clear 
pattern. This is another item to consider when building an instrument for actual use. The best 
performance can be achieved by placing a reflector or another off-axis mirror in front of the 
aperture to control the distribution and intensity of light reaching the detector. 
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Figure 40 Interference Pattern after Final Modifications to AMI 
E. TESTING WITH ULTRAVIOLET SOURCE 

The ultimate goal of this thesis was to examine the feasability of using the AMI as a high 
resolution instrument in the ultraviolet portion of the electromagnetic spectrum. The final phase 
of testing, therefore, required the use of a light source which operated in the ultraviolet 
wavelength range. After considerable thought, the mercury line at a wavelength of 2537-A was 
chosen. The mercury line offered an intense source at a wavelength close to the minimum 
wavelength acceptable for the diffraction grating of the prototype instrument. Recall from 
Equation (20) that a wavelength of 2537-A requires a grating with p min =2628 lines/mm. The 
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prototype grating had a density of 2400 lines/mm which, although not perfect, is still fairly close 
to optimum). 

Operating the AMI in the ultraviolet range presented two problems. The first problem 
was moving and adjusting the mirrors with an ultraviolet source. Since the light from the UV 
source was not visible to the naked eye, the mirrors could not be aligned visually. Because of the 
change in wavelength, the reflecting mirrors had to be moved and realigned. Initial placement of 
the mirrors was accomplished using a computer template. The template, which was taped to the 
workbench, showed the position of the incident beam, the diffraction grating, the two plane 
mirrors and the output beam to the detector for an incidence angle of 0 I =45°. 

The other problem which needed to be addressed was detecting the interference pattern. 
The computer camera did not operate in the UV range and required some modifications to detect 
light in that region of the spectrum. After careful consideration, an image intensifier was used to 
alleviate the problem. The image intensifier, a ITT model F4145 plate image intensifier, was 
bolted directly to the front of the solid state computer camera. Photons entering the intensifier 
strike a photocathode, causing electrons to be ejected and accelerated by a potential into a 
MicroChannel plate where they strike a phosphor screen. The phosphor screen emits visible light 
:h is then detected by the camera. Power to the image intensifier was provided via a separate 
circuit and was adjustable from 0 to 15 volts. 

Once modifications to the camera were complete, the mirrors were placed according to 
the template and power to the lamp and image intensifier was turned on. As expected, the zero 
order mirror was easy to align. The detector was used to check the placement of the minus one 
order mirror and showed that the template was an accurate indicator of the various reflection 
angles and mirror positions. The image intensifier seemed to work best for supply voltages 
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between 5 and 6 volts. Any power settings above 6 volts saturated the CCD camera. Final 
adjustments were accomplished by adjusting the zero and minus one order mirrors so that both 
the zero and minus one beams overlapped at the center of the detector. The setup worked so well 
that adjustments were completed for the UV source almost as fast as for the visible light sources. 
The only problem with the configuration was that the intensity of the zero order beam was five to 
ten times brighter than the minus one order beam. When power to the intensifier was adjusted for 
the minus one order beam , the zero order beam saturated the computer camera. The large 
difference in beam intensities eliminated any possibility of establishing an interference pattern. 

The incidence angle was adjusted from 45° to 5° to try and shift some of the intensity 
from the zero beam to the minus one beam and equalize the beam intensities. This new setup 
shifted the configuration of the mirrors and detector so that the minus one order mirror was on 
one side of the incident beam and the zero order mirror and detector were on the other side of the 
beam. The new configuration worked better than expected and shifted too much of the intensity 
from the zero to minus one order beam. When power was adjusted for the zero beam the 
intensity of the minus one beam saturated the camera, again preventing the appearance of an 
interference pattern. 

The next logical step was to balance the beam intensities by configuring the instrument 
for some incidence angle between 5° and 45°. The difficulty with this was that the m=l beam 
and ultraviolet emission lines in the vicinity of 3000-A (also from the mercury lamp) interfered 
with the m=0 and -1 beams at the detector for most incidence angles between 5° and 45°. As a 
result, there was no usable incidence angle for the mercury lamp that equalized the intensities of 
the beams enough to create an interference pattern. One possible solution would be to use a 
diffraction grating with a higher ruling density so that more of the intensity is concentrated in the 
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m=0 and -1 beams. Unfortunately, no other gratings were immediately available and testing of 
the AMI in the ultraviolet regi< was discontinued. 
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V. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 


A. SUMMARY OF FINDINGS 

Over the course of this thesis, numerous experiments were conducted to validate the All- 
Reflection Michelson Interferometer as a high resolution, compact instrument for use in the EUV 
and FUV portion of the electromagnetic spectrum. 

First, the instrument was examined using a mercury vapor lamp that operated at a 
wavelength of 5461-A. Test results showed that the AMI operated exactly as predicted. 
Establishing an interference pattern at the detector took a matter of minutes and the device was 
easy to tune with minor adjustments to the minus one order mirror. The resulting interference 
patterns were analyzed using IDL procedures designed to produce interference profiles and 
Fourier spectrums. As expected, the Fourier transform of the interference profiles returned a 
monochromatic spectrum. 

Next the prototype was operated using a sodium source that produced a doublet emission 
at a wavelength of 5889-A. The doublet emission was chosen to test AMI performance with a 
non-monochromatic source and determine the resolving power of the instrument. Experiments 
revealed that the instrument operated equally well with a non-monochromatic source and the 
resolution of the instrument agreed fairly closely to theoretical predictions. For reasons not yet 
understood, the resolution of the emission peaks varied depending on the number of rows used to 
sum the intensities. The greatest definition between peaks occurred when the intensities were 
summed approximately 32 to 64 rows at a time. 
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The prototype was modified for operation in the ultraviolet and was tested using a 
mercury source with a wavelength of 2537-A. Although an interference pattern was never 
established, some very useful information was gathered, including the efficacy of the image 
intensifier/CCD camera assembly for use as a detector in the ultraviolet portion of the spectrum. 

There were two reasons an interference pattern was never generated using a UV source. 
First, there was intereference from other diffraction orders and emission lines. Second, the 
difference in intensity between the beams precluded any pattern from forming. The first problem 
could be resolved by using a different source and the second problem could be eliminated by 
increasing the ruling density of the grating and determining the optimum incidence angle for the 
wavelength under observation. 

The experiments conducted during this thesis provided some extremely useful 
information and validated the use of the AMI as a high resolution sensing instrument in the 
ultraviolet region. The optical elements are simple, lightweight and relatively inexpensive. The 
instrument in general is compact, is easy to tune and adjust and very reliable. Also, since the 
AMI has no moving parts, it is sturdy and should be fairly easy to construct. 

B. RECOMMENDATIONS FOR FURTHER STUDY 

With the successful testing of a functioning prototype complete, all future work should 
be concentrated towards the development of an actual flight instrument. The engineering 
development should focus on a number of important elements. First, different methods of 
collecting and collimating the incident light should be examined and tested. Second, the type and 
optimum grating density of the diffraction grating should be explored in more detail. Third, the 
methodology for mounting the optical elements to reduce vibrations due to shock and 
temperature fluctuations during flight should be investigated. Fourth, detectors suitable for 
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sensing interference patterns in the EUV and FUV region of the electromagnetic spectrum 
should be examined. Finally, a suitable arrangement (including dimensions, weight and 
placement of elements) to conform to sounding rocket or space vehicle payload constraints 
should be explored. 
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APPENDIX A 


IDL PROCEDURES 


PRO INITVAL,T,WLSIG,WLOFF,WLCEN,IO,XSEC 

; THIS PROCEDURE GIVES INITIAL VALUES FOR : HE 
; OXYGEN 1304 LINE PROFILE PROCEDURES 


WLSIG=.0478D 

WLOFF=.0478D 

PEAKS(ANGSTROMS) 

WLCEN=1304.86D 

I0=1D 

A=15.9994D 

F12=.0485D 

XSEC=1.16D-14*WLCEN*F 12*((A/T) A .5) 

RETURN 

END 


; LINE DISPERSION(ANGSTROMS) 

; OFFSET OF EMISSION 

CENTER WAVELENGTH(ANGSTROMS) 
; AMPLITUDE 

; ATOMIC WEIGHT OF OXYGEN 
; ATOM(GRAMS) 

; OSCILLATOR STRENGTH 
; RESONANT SCATTERING CROSS 
; SECTION(CM A 2) 


PRO SOLARLIN,T,SOLENT, WL 

; THIS PROCEDURE CALCULATES THE SOLAR OXYGEN 1304.86 LINE PROFILE 

INITV AL,T,WLSIG,WLOFF,WLCEN,10,XSEC 

NPTS=2000 

WL=DBLARR(4000) 

SOLINT=DB LARR(4000) 

WL=WLCEN -. 2 5+FINDGEN (2 * NPTS)/ (NPTS * 4) 

SOLINT=IO*(EXP(-((WL-WLCEN+WLOFF) A 2)/WLS1G A 2)+EXP(-((WL-WLCEN- 
W LOF F) A 2)/WLSI G A 2)) 

RETURN 

END 


PRO TERLIN,WL,TERESINT 
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; THIS PROCEDURE CALCULATES THE TERRESTRIAL 1304.86 LINE USING A 
; DOPPLER LINE PROFILE 

NPTS=2000 
WLCEN=1304.86D 

WL-WLCEN-.25+FINDGEN(2*NPTS)/(NPTS*4) 

WLSIG=1.862E-6*WLCEN 
TERESINT=EXP(-((WL-WLCEN) A 2)/WLSIG A 2) 

RETLTRN 

END 


PRO VOITLIN,T,WL,VOIGTINT 

; THIS PROCEDURE CALCULATES THE TERRESTRIAL 1304.86 LINE PROFILE BASED 
; ON THE VOIGT FUNCTION 


A=1.108D-2 

WLCEN=1304.86D 

FREQ=DBLARR(4000) 

U=DBLARR(4000) 

H=DBLARR(4000) 
VOIGTINT=DBLARR(4000) 
DOPWlDTH=2.47D8*SQRT(T) 
FREQCEN=2.998D8/(WLCEN* ID-10) 
FREQ=2.998D8/(WL* 1D-10) 
U=(FREQ-FREQCEN)/DOPWIDTH 
H=VOIGT(A,ABS(U)) 
VOIGTINT=H/(DOPWIDTH A .5)*(!DPI A .5) 

RETURN 

END 


PRO VOBSINT,T,COLDEN,WL,VOBSINT 
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; THIS PROCEDURE CALCULATES THE OBSERVED OXYGEN 1304.86 LINE PROFILE 
; USING THE VOIGT FUNCTION TO MODEL THE TERRESTRIAL LINE 

INITVAL,T,WLS1G,WLOFF,WLCEN,IO,XSEC 
SOLARLI1 ,T, WLSIG. WLOFF.WLCENJO.SOLINT, WL 
VOITLIN,T,WL,VOIGTINT 
VOBSINT=DBLARR(4000) 

VOBSINT=SOLINT*EXP(-(COLDEN*XSEC*VOIGTINT)) 

RETURN 

END 


PRO INTVSCOL,T,COLDEN,INT 

; THIS PROCEDURE CALCULATES THE NORMALIZED INTENSITY OF THE 
; OBSERVED OXYGEN 1304.86 LINE VERSUS AN ARRAY OF COLDENS 

INITVAL,T,WLSIG,WLOFF,WLCEN,IO,XSEC 
COLDEN=FLTARR( 101) 

INT=FLTARR(101) 

COLDEN(0)= 1E16 

FOR J=1,100 DO BEGIN 
COLDEN(Jf)=COLDEN(J-1 )* 1.26 
ENDFOR 

FOR J=0,100 DO BEGIN 
VOBSINT,T,COLDEN(J),WL,VOBSINT 
INT(J)=TOTAL(VOBSINT) 

ENDFOR 

INT=INT/INT(0) 

RETURN 

END 


PRO INTVCOL1 
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; THIS PROCEDURE PLOTS THE GROWTH CURVE OF THE OXYGEN 1304.86 LINE 
; SEVERAL TEMPS AND SENDS THE DATA TO A CGM FILE 

SET_PLOT,'CGM' 

DEVICE,FILENAME='GROWTH.CGM' 

INTVSCOL,300,COLDEN,INT 

PLOT_OI,COLDEN,INT,XSTYLE=l,XRANGE=[lE16,lE24],$ 

XT1TLE-COLUMN DENSITY (CM A -2)',$ 

YTITLE='INTENSITY RATIO (lobserved/Isolar)' 

OPLOT,[lE17,lE18],[.35,.35] 

XYOUTS,5E18,.35,'T=300K',SIZE=75 
OPLOT, [ 1E17,1E18], [,3, .3 ] ,LINE= 1 
XYOUTS,5E18,.3,'T=600K',SIZE=.75 
0PL0T,[IE17,1E18],[.25,.25],LINE=2 
XYOUTS,5E18,.25,'T=900K',SIZE=.75 
OPLOT,[1E 17.1E18],[.2,.2],LINE=3 
XYOUTS,5E18,.2,'T=1200K',SIZE=.75 
OPLOT,[ 1 El 7,1 El 8],[. 15,. 15],LINE=4 
XYOUTS.5E 18,.l 5,'T=1500K',SIZE=.75 
OPLOT,[1E17.1E18],[.l,.l],LINE=5 
XYOUTS,5E 18,. 1 ,'T= 1800K',SIZE=.75 

FOR T=600,l 800,300 DO BEGIN 
INTVSCOL,T,COLDEN,INT 
OPLOT,COLDEN,INT,LINE=(T-300)/300 
ENDFOR 

DEVICE,/CLOSE 

SET_PLOT,'WIN' 

RETURN 

END 


PRO IMAGE,FILE,IMAGE 



; THIS PROCEDURE READS TIFF FILES PRODUCED BY THE EDC CAMERA AND 
; SENDS THE PICTURE TO THE SCREEN 

FILENAME='C:\EDC\'+FILE 

IMAGE=TIFF_READ(FILENAME) 

LOADCT.6 

TV,IMAGE,4,8 

RETURN 

END 


PRO PROFILE,FILE,PROF 

; THIS PROCEDURE READS TIFF FILES PRODUCED BY THE EDC CAMERA AND 
; SUMS THE INTENSITIES BY COLUMN TO PRODUCE AN INTERFERENCE PROFILE 

FILENAME-C :\EDCY+FILE 
IMAGE=TIFF_READ(FILENAME) 

PROF=TOT AL(IMAGE,2) 

PROF=PROF-TOTAL(PROF* 1 .)/753. 

RETURN 

END 


PRO SPECTRUM,FILE,PROF,SPEC 

; THIS PROCEDURE READS TIFF FILES PRODUCED BY THE EDC CAMERA, SUMS 
; THE INTENSITIES BY COLUMN TO PRODUCE AN INTERFERENCE PROFILE AND 
; TAKES THE FOURIER TRANSFORM 

FILENAME-C:\EDCY+FILE 
IMAGE=TIFF_READ(FILENAME) 

PROF=TOTAL(IMAGE,2) 

PROF=PROF-TOTAL(PROF* 1 ,)/753. 

SPEC=ABS(FFT(PROF,-1)) 

RETURN 

END 


PRO GETSPEC,FILENUM,SPEC 



; THIS PROCEDURE READS DATA FROM A TIFF FILE USING THE IMAGE 
; PROCEDURE, SUMS THE DATA NUM ROWS AT A TIME, TAKES THE FOURIER 
; TRANSFORM OF THE PROFILE AND STORES THE DATA AS A RECORD IN THE 
; FILE NEWSPEC.DAT 

SPEC=F LT ARR( 753) 

FDAT=FLTARR(753) 

IMAGE,FILE,IMAGE 
FINISH=F1X(244./NUM)-1 

FOR J=0,FINISH DO BEGIN 
DAT=TOTAL(IMAG(*,NUM*J:NUM*J+(NUM-1 )),2) 

FD AT=ABS(FFT(D AT-TOTAL(DAT)/753.,-1)) 

SPEC=SPEC+FDAT 

ENDFOR 

IF NUM*J+(NUM-1) LE 244 THEN BEGIN 
DAT=TOTAL(IMAG(*,NUM*J:*),2) 

FD AT=ABS(FFT(DAT-TOTAL(DAT)/7531)) 

SPEC=SPEC+FDAT 

ENDIF 

RETURN 

END 
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APPENDIX B 


EDC FILES USED IN FIGURES 


1MAGE8.TIF 

IMAGE12.TIF 

IMAGE14.TIF 

IMAGE15.TIF 

IMAGE26.TIF 

IMAGE32.TIF 

IMAGE39.TIF 


Figure 24, Figure 25, Figure 26 
Figure 30 

Figure 27, Figure 28, Figure 29, Figure 30 
Figure 30 

Figure 31, Figure 32, Figure 33 
Figure 34, Figure 35, Figure 36, Figure 37 
Figure 40 
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